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Abstract. In the classical vacuum Maxwell-Lorentz theory the self-force of a charged 
point particle is infinite. This makes classical mass renormalization necessary and, 
in the special relativistic domain, leads to the Abraham-Lorentz-Dirac equation of 
motion possessing unphysical run-away and pre-acceleration solutions. In this paper we 
investigate whether the higher-order modification of classical vacuum electrodynamics 
suggested by Bopp, Lande, Thomas and Podolsky in the 1940s, can provide a solution 
to this problem. Since the theory is linear, Green-function techniques enable one 
to write the field of a charged point particle on Minkowski spacetime as an integral 
over the particle’s history. By introducing the notion of timelike worldlines that are 
“bounded away from the backward light-cone” we are able to prescribe criteria for the 
convergence of such integrals. We also exhibit a timelike worldline yielding singular 
fields on a lightlike hyperplane in spacetime. In this case the held is mildly singular 
at the event where the particle crosses the hyperplane. Even in the case when the 
BoppPodolsky held is bounded, it exhibits a directional discontinuity as one approaches 
the point particle. We describe a procedure for assigning a value to the held on 
the particle worldline which enables one to dehne a hnite Lorentz self-force. This 
is explicitly derived leading to an integro-differential equation for the motion of the 
particle in an external electromagnetic held. We conclude that any worldline solutions 
to this equation belonging to the categories discussed in the paper have continuous 
4-velocities. 


Keywords'. Self-force, Radiation reaction, Higher-order electrodynamics, Bopp-Podolsky 
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1. Introduction 

For many applications it is reasonable to model moving charges in terms of classical 
charged point particles. In accelerator physics, for example, it is usually neither desirable 
nor feasible to model particle beams in terms of extended classical charged bodies or 
of quantum matter. Therefore a mathematically consistent theory of classical charged 
point particles is of high relevance. 
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Unfortunately, such a theory does not exist so far. Of course, there is no problem 
as long as we restrict to a classical charged test particle and neglect the particle’s 
self-interaction. Then the equation of motion is just the relativistic generalization of 
Newton’s equation of motion with the rate of change of particle momentum equated 
to the Minkowski (relativistic) Lorentz force in a given external field and everything is 
fine. If, however, the self-field is taken into account, the theory becomes pathological. 
According to the Maxwell-Lorentz theory in vacuo, the electromagnetic field of a point 
charge becomes infinite at the position of the charge, so the particle experiences an 
infinite self-force. This infinity is so bad that the field energy in an arbitrarily small 
ball around the source is infinite which leads to an infinite term in the equation of 
motion of a point charge. Dirac (1938) suggested to counter-balance this infinity by 
postulating that the point charge carries a negative infinite “bare mass” which leads to 
the Abraham-Lorentz-Dirac equation. Even if one accepts this ad-hoc idea of an infinite 
bare mass, the problem has not been solved. The Abraham-Lorentz-Dirac equation is 
known to possess unphysical behavior such as run-away solutions and pre-acceleration. 
Reviews of this dilemma, including detailed accounts of the history, can be found in the 
comprehensive monographs by Rohrlich (2007) and Spohn (2007). 

Since the Maxwell-Lorentz theory with point charges in vacuo does not lead to 
a consistent equation of motion of charged point particles, one might think about 
modifying this theory. In the course of history at least two such modifications have been 
suggested, both motivated by the desire of solving the problem of an infinite self-force, 
namely the Born-Infeld theory and the Bopp-Podolsky theory. The Born-Infeld theory 
is by far the better known of the two. This theory, which was suggested by Born & Infold 
(1934), modifies the source free Maxwell vacuum theory by introducing non-linearities 
containing a new hypothetical constant of Nature b with the dimension of a (magnetic) 
field strength. For b —> oo the Maxwell equations in vacuo are recovered. The fact 
that the Maxwell theory is very well verified by many experiments is in agreement with 
the Born-Infeld theory as long as b is sufficiently large. By contrast, the Bopp-Podolsky 
theory retains linearity but introduces higher-derivative terms proportional to a factor £ 2 
where £ is a new hypothetical constant of Nature with the dimension of a length. Again, 
for £ —y 0 the Maxwell-Lorentz theory is recovered. The Bopp-Podolsky theory was first 
suggested by Bopp (1940). It was independently rediscovered by Podolsky (1942). Both 
Bopp and Podolsky formulated their theory in terms of an action functional and then 
derived the field equation which is of fourth order in the electromagnetic potential. As 
noted by both Bopp and Podolsky, this fourth-order equation is equivalent to a pair of 
second-order equations in a certain gauge. If rewritten in this form, the Bopp-Podolsky 
field system coincides with those of a theory suggested by Lande & Thomas (1941). 
Similar to the Born-Infeld theory, the Bopp-Podolsky theory was first formulated as 
a classical field theory but with the intention of deriving a quantum version later. In 
particular, Podolsky pursued both the classical and the quantum aspects of the theory 
in several follow-up articles with different co-authors, see Podolsky & Kikuchi (1944), 
Podolsky & Kikuchi (1945) and Podolsky & Schwed (1948). In the present article we 
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are interested only in the classical theory. 

In both the Born-Infeld theory and the Bopp-Podolsky theory the self-held is 
bounded for a static point charge, i.e., for a point charge that is at rest in some inertial 
system in Minkowski spacetime. This was shown already in the earliest articles on these 
theories. Moreover, in both theories for such a charge the held energy in a ball of radius 
R around the charge is finite, even in the limit R —> oo. To the best of our knowledge 
in the Born-Infeld theory little is known regarding such finiteness for accelerated point 
charges. In the Bopp-Podolsky theory the only result about accelerated point charges 
that we are aware of is due to Zayats (2014), who showed that the self-force is finite for 
a uniformly accelerated particle on Minkowski spacetime. 

It is the purpose of this paper to add some results on the finiteness of the self-force 
in the Bopp-Podolsky theory. In our view, these results give strong support to the 
idea that the Bopp-Podolsky theory provides a consistent theory of classical charged 
point particles including the self-force. In Section [2] we briefly review the basic field 
equations of the Bopp-Podolsky theory on Minkowski spacetime, emphasizing the fact 
that because of their linearity Green-function techniques can be used. In Section [3] we 
restrict the equations of Section [2] to the case where the source of the electromagnetic 
field is a point charge with a prescribed worldline on Minkowski spacetime. We discuss 
various ways of writing the field strength at a point off the worldline as an integral over 
the particle’s history. As a first example, we treat the simple case of a point charge that 
is at rest in an inertial system. In Section [4] the general problem of assigning a value of 
the field strength on the worldline is discussed. This is a precursor to the formulation of 
an integro-differential equation for the motion of a point charge that accommodates its 
finite self-force in an external electromagnetic field. As a second example, we treat the 
case of a uniformly accelerated charge. In Section [5] we present our main results on the 
finiteness of the field and of the self-force for general motion. We show that the self-force 
is finite unless the worldline approaches the light-cone in the past in a very contrived 
manner. As a third example, we discuss such a pathological worldline where the self¬ 
field actually diverges on a lightlike hyperplane and the self-force diverges at one point 
on the worldline. However, we demonstrate that even in this case the singularity of the 
field is so mild that it does not cause a problem for the equation of motion. In Section [6] 
we discuss how the Abraham-Lorentz-Dirac equation comes about in a particular limit 
as i —)■ 0, after classical mass renormalization. 

In the body of the paper we formulate the Bopp-Podolsky theory on Minkowski 
spacetime in an inertial coordinate system. However, we have added an appendix 
where we consider the Bopp-Podolsky theory on a curved spacetime in arbitrary 
coordinates. This allows us to derive the dynamical (Hilbert) electromagnetic stress- 
energy-momentum tensor of the theory using exterior calculus. The appendix also 
includes a derivation of the relativistic Lorentz force in the Bopp-Podolsky theory from 
this tensor which is crucial for our reasoning in the body of the paper. 
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2. Bopp-Podolsky theory 

We consider a time and space oriented Minkowski spacetime with standard inertial 
coordinates x = (x°, x 1 , x 2 , x 3 ), with metric tensor 

g = Vabdx a ®dx b (1) 

where (r] a b) = diag(—1,1,1,1). In this article all tensor field components on Minkowski 
spacetime are with respect to the class of global parallel bases adapted to these 
coordinates. Members of this class are related by elements of the proper Lorentz group, 
SO(3, 1). Here and in the following, Einstein’s summation convention is used for latin 
indices which take values 0,1, 2, 3 and for greek indices which take values 1, 2, 3. Latin 
indices are lowered and raised with r) ab and with its inverse rf b , respectively. We use 
units in which the the speed of light c equal to 1. 

The higher-order electrodynamics suggested by Bopp (1940) and, independently, 
by Podolsky (1942) is based on the gauge invariant action functional 

S[A] = {-^ Fab Fab + ^-d c F ab d c F ab - A a f^j d A x (2) 

where A4 is some compact region of Minkowski spacetime yielding a finite S'[A]. Here 
A a is the electromagnetic potential 

F ab = d a A b - d b A a (3) 

is the electromagnetic field strength, j a is a conserved current density 4-vector field, 
{d a = d/dx a } and i is a hypothetical new constant of nature with the dimension of 
a length. It proves expedient to derive the Bopp-Podolsky field equations in terms 
of smooth fields and a smooth current source on A4 and then discuss particular 
singular solutions associated with a source having support on a tinrclike worldline. This 
eliminates the need to perform variations of ([2]) in a distributional context. Note that, 
for deriving the field equations by variational methods, the action functional (J2J) can be 
equivalently replaced with 

S[A] = f M (l h FabFah + ^ daF^Fcb - A 0 f^J d A x (4) 

because the integrands differ only by a total divergence. 

The field equations of the Bopp-Podolsky theory result from varying the action 
functional (]2]) or (J4]) with respect to the potential. They read 

d b F ba - fnd b F ba = -4vr j a (5) 

or, in terms of the potential in the Lorenz gauge d b A b = 0, 

UA a -i 2 U 2 A a = -47 Tj a (6) 

where □ = d b d b is the wave operator. Field equations involving the operator D 2 have 
also been investigated by Pais & Uhlenbeck (1950), cf. Pavlopoulos (1967). 
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Both Bopp and Podolsky observed that the fourth-order differential equation flhl) 
can be reduced to a pair of second-order differential equations. More precisely, (J6]) is 
equivalent to 

□4* = — Aitj a , d a A a = 0, (7) 

□i a - r 2 A a = - A it j a , d a A a = 0. (8) 

This can be demonstrated in the following way. Assume we have a solution A n to (161) . 
Then we define 

A a = A a - £ 2 DA a , A a = -fDA a (9) 

and it is readily verified that (JTJ) and OH]) are indeed true. Conversely, assume that we 
have solutions A a and A a to (J7]) and OH]), respectively. Then we define 

A a = A a - A a (10) 

and it is readily verified that (ED is true. This gives a one-to-one relation between 
solutions to (E]) and pairs of solutions to (JZ]) and OH]) which allows one to view the Bopp- 
Podolsky theory as equivalent to a theory based on the two equations (ED and OH]) • The 
latter was suggested, shortly after Bopp but independently of him and shortly before 
Podolsky, by Lande & Thomas (1941). In a quantized version of the Lande-Thomas 
theory, (ED describes the usual (massless) photon while (JH]) describes a hypothetical 
“massive photon” whose Compton wave length is equal to the new constant of nature 


To find the retarded solution to the fourth-order Bopp-Podolsky equations (ED, for 
any given divergence-free source j a , we use the reduction to the second-order equations 
(ED and © and then apply standard Green-function techniques. The details were worked 
out by Lande and Thomas, see Section 9 in (Lande & Thomas 1941). The retarded 
solution to 

DG{x - y) = -8(x - y) (11) 

is 

( (2vr )~ 1 6(D(x - y) 2 ) if y < x , 

G(x -y) = l 

y 0 otherwise, 

and the retarded solution to 

□G (x -y) - t~ 2 G(x -y) = -S(x-y) 


( 12 ) 

(13) 


is 


G(x-y) 


(2tt) l 5(D(x — y) 2 ) 


Ji(£-~ 1 D(x — y)) 

A%£D(x-y) V X ' (14) 


0 


otherwise. 
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Here and in the following, J n is the nth order Bessel function of the first kind, y < x 
means that y is in the chronological past of x and D(x — y) is the Lorentzian distance 
between these two events, 

D(x-y) = \/-r]ab{x a - y a )(x b - y b ) . (15) 


Hence, the retarded solution to (J7j) is 

A a (x) = 4 tt f G(x — y) j a (y) d A y , 

J R 4 

the retarded solution to (jS]) is 

A a (x) = 4vr [ G(x - y) j a (y) d l y , 

J R 4 

and the retarded solution to (JH]) is 

A a (x) = 4v rf (d(x - y) - G(x - y)^j j a (y) d 4 y 

■3a(y) <&V- 


(16) 


(17) 


J^- l D(x-y)) . 


(18) 


Jy<x tD(x - y) 

The Lorenz gauge condition is satisfied by A a , A a and A a when the current density 
satisfies the continuity equation 


daf = 0. (19) 

Mathematically A a is the retarded potential of the standard Maxwell theory. For i —y 0 
we have G —> 0 and hence A a —y 0, so in this limit the standard Maxwell theory is 
recovered, as is obvious from (|2il . 

Equation (fT8l) can be viewed as a map that assigns to each current density j 
the corresponding retarded Bopp-Podolsky potential A. A general framework for 
investigating the question of whether this map is well-defined would be to assume that 
j is a (tempered) distribution and to ask if A is again a (tempered) distribution. As 
the Green function G — G does not satisfy a fall-off condition in all spacetime directions, 
this is a non-trivial question. In this paper we restrict to a more specific question. We 
assume that j is the current density associated with a point charge, and investigate 
whether the integral on the right-hand side of (fTHj) converges for events x in the future 
of the worldline of the charge. If this is the case, the left-hand side of (1T81) is, of course, 
well defined, not only as a distribution but even as a function. 

Up to now we have discussed only how the electromagnetic field can be calculated 
from its source, i.e., from the current that generates the field. We also need the 
equation for the Minkowski Lorentz force density which arises from the divergence 
of an electromagnetic stress-energy-momentum tensor T a b. In inertial coordinates on 
Minkowski spacetime, it reads 

d c T ca = F abJ b . (20) 


This equation defines the Minkowski force density that the field F a b exerts on the current 
j b . In the case that the current is concentrated on a worldline, it gives the self-force 
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Figure 1. When the future pointing timelike worldlinc £(r) approaches a lightlike line 
(shown dashed) asr-> r m i n , its future is the subset of Minkowski spacetime, partially 
shown in grey above, bounded by a lightlike hyperplane containing the dashed line. 
The grey domain only consists of events which can be reached from the worldline along 
a future-oriented timelike curve. 

in terms of field components on the worldline. A derivation of (l20j) in terms of the 
stress-energy-momentum tensor T a f, associated with the Bopp-Podolsky theory, is given 
in the Appendix. 

3. The field of a point charge off the worldline 

Let £(r) = (^°(r),^ 1 (r),^ 2 (r),^ 3 (r)) be an inextendible timelike C°° curve parametrized 
by proper time, 

mt(T) = -i. ( 21 ) 

Inextendible curves have parametrizations with r belonging to an open interval 
] Tmin, r max [ where r min elU {-oo} and r max e R U {+oo}H 

Consider the future of £(r), i.e., the set of all events that can be reached from the 
worldline along a future-oriented timelike curve. If the worldline approaches a light-cone 
asymptotically for r —> T min , its future is bounded by a lightlike hyperplane, see Fig. HI 
otherwise it is all of R 4 . To each x, in the future of £(r), but not on the worldline, we 
assign the retarded time r R (x), defined by the properties that 

(x a - C(t r (x))^J (x a - £ a (r R (x))) = 0 , x° > ^°(t r (x)) , (22) 

and introduce the retarded distance 

r R (x ) = -C(r R (x))(x a -^a(r R (x))^ . (23) 

Unless otherwise specified, in the following, all events x will belong to the future of the 
worldline. However, whether such events lie on or off the worldline will be made explicit. 
For events x off the worldline one has 

x a = £ a (r R (x)) + r R (x)(i a (r R (x)) +n a (x)^j (24) 

Some curves with unbounded acceleration reach past or future infinity in a finite proper time. 


with a well-defined spatial unit vector n(x), 

n a {x)n a (x) = 1, £ a (r R (x))n a (x) = 0, 

see Fig. [2j By differentiation, we find 

d b T R {x ) = -6(r R (a?)) - n b {x) 

and 

d b r R (x) = n b (x ) + r R (x)£ a (T R (x))n a (x) (i b {r R (x)) + n b (x) ) . (27) 


(25) 

(26) 


For a sequence of events x N that approaches the worldline, x N —> £(r 0 ) as N —> oo, 
t r (xn) — > To and t r (xn) — > 0. The limits for n(x at), d b r R (x n) and d b r R (x jv) do not 
exist as N —> oo. 

We model a point charge with worldline £(r) by the distributional current density 


Ja{x) = q 


&( x - £M) L(t) dr 


(28) 


where q is its electric charge. Then (fT6l) gives the standard Lienard-Wiechert potential, 

qL(r R (x)) 


AJx) = 


r r( x ) 


and (jT8|l reads 


AAx] = 


rr R ( X) J^i-'Dix - Zir))) 


L(t) dr 


(29) 


(30) 


t„,„ D(x - Z(t)) 

Both (|29l) and (l30l) are defined for events x off the worldline. In Section [4] we investigate 
what happens to the potentials and fields as x approaches the worldline. 

Note that, in contrast to the Lienard-Wiechert potential (|29l) . the potential (|30l) at 
an event x depends on the whole history of the charge from r = r m ; n up to r = t r (x). 
(The same is true, in general, in the standard Maxwell-Lorentz theory with point charges 
on curved spacetimes.) The integrand in (1301) is bounded for r —> t r (x) because 


j.^.-ew)) jmx )) 


U(r) 


for r —> t r (x) 


(31) 


D(x-£(t)) ' 2« 

where we have used the Bernoulli-l’Hopital rule and T((0) = 1/2. By contrast, for 
t — > T m ; n , the individual components ^ a (r) may blow up arbitrarily. Therefore, the 
existence of the integral on the right-hand side of (l30l) is not guaranteed. We show later 
that, for a fairly large class of worldlines, this integral does converge even absolutely, as 
a Lebesgue integral or as an improper Ricmann integral, for all x in the future of the 
worldline; however, we also give a (contrived) example where it does not converge for 
some x in the future of the worldline. 

On the assumption that the integral converges, at all events x we can differentiate 
()30l) off the worldline with respect to x b . Antisymmetrizing the resulting expression 
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gives the field strength 


q 


Fab(x) = [£b(T R (x))n a (x) - £ a (T R (x))n b (x] 


(32) 


q_ 

i 2 


J 2 (l l D(x - £(r))j 


n , x {xb^ib{r))i a {T)-{x a -ia{r))i b {T))dT. 

D{x~£{t)) 

Here we have used (125]) and the identities 2 J\{z) = z(J 0 (z) + J 2 (^)) and 2 J[(z) = 
Jo(z) — J 2 (z) of the Bessel functions. Again, we postpone the discussion of what happens 
if x approaches the worldline to Section [4] 

We observe that, keeping x fixed, we may use ( = D(x — £(r)) as the parameter 
along the worldline. Indeed, differentiation of the equation 

C 2 = -(x a -f a (r))(x a -f a (r)) (33) 

yields 

(d( = i a (r)(x a ~ ia{T))dr. (34) 

As, by the reverse Schwartz inequality for timelike future-oriented vectors, 

i a ir){x a - f a (r)) < -C < 0 for 0 < ( < oo , (35) 


is monotonically decreasing along the worldline. This guarantees that the equation 
( = D(x — £(r)) can be solved for r, 

C = £>(&-£(?)) 


T = W(C,X) . 


(36) 
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As proper time r runs from r min to t r (x), the new parameter ( runs (backwards) from 
oo to 0. Hence at all events x (130]) can be rewritten as 


A (x) = - q r wmut) 

flU Wo i\r){x b -Ur)) 


d(. 


(37) 


T=Vj(C > ,x) 

Note that, if we regard x as a, parameter, (l37l) has the form of a Hankel transform which 
transforms a function of £ to a function of \ji. 

We may use the parameter ( in the formula for the field strength as well. If such 
a change of the integration variable is performed on the right-hand side of (]32|) . the 
resulting equation reads 


F ab(x) = ^ (6(rB(*))no(*) - £a(TR(x))n b (x) 


°° ( (x b - ^b(r))ia(r) - (x a - £a(r))£b(T) 


(38) 


Q_ _ 

^ Jo £ c (r)(x c — £ c (r)) 


C 


T=Zj(C,,x) 

An alternative expression for the field strength is obtained if, off the worldline, we 
differentiate (157|) with respect to x b and antisymmetrize, 


q r°° J l((/t)[£a(T){x b -£b(T)) - €b(T)(x a - £a(T)) 

F ab( x ) = V 

t - /n 'Z c (t)(x c -Ut)) 


d( 


T=Zj(£,x) 

r-oo Ji{C/()(ia{r)(x b - ib{r))-i b {r)(x a - C a (r))) (l + i d {T)(x b - €b(r)) 




£ c (t)(x c ~ £ c (t 


(39) 


d(- 


T=Zu((^,x) 


can be derived directly from (|38|) by integrating its second term by 


Alternatively 
parts. 

It is worth noting from (1371) that we can derive another form of the potential by 
performing an integration by parts and using the identity — J\ = J' Q of Bessel functions. 
The resulting equation 

q£a{r R (x)) 


A a (x) = 


tr(x) 


(40) 


£a(r) - 


- q 


£a{x) (^1 + £ d (r)(x d — £ d (r 

£ C (t)( x c ~ £c( r )) 


£ b (r)(x b — £b(r)) 


M(/t)(d( 


T=Vj(C > ,x) 


gives the deviation of the potential from the Lienard-Wiechert potential. With the 
help of standard asymptotic formulas for the Bessel function Jq the right-hand side can 
be written as a power series in l. Such asymptotic (i.e., in general non-convergent) 
expansions have been used, e.g., by Frenkel (1996) (also see Frenkel & Santos (1999)) 
and Zayats (2014). 
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Example 1: Charge at rest 


The simplest case one can consider is a charge at rest in an appropriately chosen inertial 
system; this is equivalent to saying that the worldline of the charge is a straight timelike 
line, 


£°(r) =V a r (41) 

with a constant four-vector V satisfying V a V a = —1. In this case (frill is a quadratic 
equation, 

( 2 = t 2 - 2 r(r R (x) - t r ( x)) + 2 r R (x)r R (x) + t r (x ) 2 , (42) 


which can be easily solved for r, 

uj((,x) = r R (x) + r R (x) - VC 2 + r R (x) 2 . (43) 


Then (13711 simplifies to 

A a (x) = - 

which is finite for all x, 


qV a f°° _ qV a 

^ Jo \/C 2 + r R (x ) 2 r R (x) 


(l - e ~ r ^ x)/ ^ (44) 


AJx) 


qV a 


for r R (x) — y 0 . 


(45) 


By differentiation of (144)1 off the worldline, or equivalently by evaluation of (138)1 or (1391) 
off the worldline, we get the held strength 

F ab (x) = V^xf { Ua ^ Vb ~ n ^ x ) V J) ( 1 - e ~ rR{x),i ~ e ~ rRix)/e ) 

= ^2 - M x )Va ) ( 1 + 0(rJi(x))) . (46) 

This expression is not defined on the worldline; if a point on the worldline is approached, 
the limit of some components depend on the direction and we say that the held displays a 
directional singularity. Although for events x on the worldline the held Q38)) is undehned, 
the held ()39l) can be evaluated there and yields F a b(x ) = 0 for all a,b. This value also 
arises by a certain averaging procedure described in the next section. 

In the rest system of the charge we have — 0 for /a — 1, 2, 3 and r R (x) = 
\J (x 1 ) 2 + ( x 2 ) 2 + (x 3 ) 2 is just the ordinary radius coordinate. In this coordinate system 
(146|) gives a radial electrostatic held with modulus 



In contrast to the Coulomb held of the standard Maxwell theory, the Bopp-Podolsky 
E(r ) stays finite as the worldline is approached. This result played a crucial role in 
the original work of Bopp (1940) and Podolsky (1942). It has the consequence that, 
at least for a charge at rest, the total held energ}@ is finite. Note, however, that the 
electric vector field cannot be continuously extended into the origin, because of the 
above-mentioned directional singularity, see Fig. |3j 


+ The total field energy £ = f R3 Tgod 3 x where T a b is defined in the appendix. 
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E x (x) 



Figure 3. The ^-component of the electric field strength on the x-axis of a charge at 
rest in the Bopp-Podolsky theory (solid) and in the standard Maxwell theory (dashed) 

4. Field of a point charge on the worldline and self-force 

The two expressions (1301) and (1371) for the potential are equivalent off the worldline, 
and so are the two expressions (138]1 and (l39|l for the field strength. We now discuss 
their behavior on and near the worldline. This is crucial because the value of the 
field strength on the worldline will determine the self-force. Both representations of the 
potential (j3H and (j37l) are well defined and continuous on the worldline. If the potential 
were differentiable, its derivative would give the field strength on the worldline without 
any ambiguity. However, it is not differentiable on the worldline. This is the reason 
why the expression (J38]), which results from differentiating (J30]), and the expression (j39lh 
which results from differentiating (1371) . behave differently on the worldline. 

First observe that (T38l) is not defined at events on the worldline because it involves 
the derivative O^tr^x) which is not defined at such events. By contrast, (1391) does 
not involve this derivative and gives a unique value for the field strength on the 
worldline, provided that the integrals converge. However, this value does not result 
from differentiating the potential on the worldline because in order to differentiate (1371) 
one must pass the derivative under the integral; this is only valid if the integrand has a 
continuous derivative, which does not occur if x is on the worldline. 

As noted, the field F a b(x ) given by (1381) or (1391) is not continuous in a neighborhood 
of events on the worldline: if a point x on the worldline is approached some of the 
components have non-zero finite limits that depend on the direction taken. This 
is manifest in (1381) where the discontinuity arises solely from the first term. Such 
behavior contrasts with similar problems in defining the self-force on a charged particle 
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in Maxwell-Lorentz electrodynamics where a similar directional dependence arises as 
one approaches the particle worldline. However, in that theory the same limiting values 
are infinite. 

Since the self-force (to be defined below) depends on the values of F ab (x ) on the 
worldline, one must assign them specific values. Such assignments should be based on 
physical criteria beyond the mathematical analysis thus far. For example an isolated free 
point charge should remain in “inertial” motion in the absence of external fields, i.e. the 
self-force should be zero. From example 1, this can be achieved by assigning the value 
zero to all F ab (x) on the inertial worldline in this case. This assignment is equivalent 
to either adopting (]39l) directly for F ab (x ) or applying a “directional-averaging” (see 
below) to (l38j) . Although for arbitrary motion one may always define a directional- 
averaging such that the value given by (]39l) arises from such a process applied to (]38l) . 
this is by no means a unique procedure. In general one may consider more complex 
averaging procedures involving constructions based on extensions of the natural Frenet 
frame defined by the worldline. 

A procedure based solely on the worldline’s tangent vector, rather than the full 
geometry of its Frenet frame can be done most easily if we introduce the (spatial) 
retarded 2-sphere, see Fig. [U 

S(T 0 ,r 0 ) = {x e M ’l t r (x) = r 0 , r R {x) = r 0 } . (48) 

The directional-average of any spacetime tensor’s inertial components K a '" b ...( x ) over 
the 2-sphere S(r 0 ,r 0 ) is defined as 

[ [ K a - b ..{x) dS (49) 

47r J Js(T 0 ,r 0 ) 

where dS is the natural surface measure on S(tq, r 0 ), induced by the ambient Minkowski 
metric. This integral depends only on To, r 0 , £(ro) and £(ro) and is manifestly Lorentz 
covariant with respect to global Lorentz transformations. I.e. if 

K a '- v Jx) = A a ' a • • • A/ • • • K a - b Jx) (50) 

then 

K“'= A a ' a ---A b , b ---K^Z. (51) 
where A a ' a <E SO( 3,1). 

Coordinating S(jo, r 0 ) with spherical polars the measure dS = sin id did dp 

and (jj9j) becomes 

2 r 27T P7T 

K a "' b ... = — / / K a "' b „Xx) smtidfidip (52) 

J o J 0 

In terms of an orthonormal tetrad (eo, ei, e. 2 , 63 ) with eo = £(ro) 

n a (x ) = cos (p sin id + sin sin d e ? 2 + cos'd (53) 

for x G S(r 0 ,r 0 ). Hence n a = 0. Since £, a (T R (x)) = ^“(to) is constant on S'(to, r 0 ) the 
hrst term on the right-hand side of (155]) averages to zero. Hence its limit for r Q —y 0 is 
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Figure 4. This figure illustrates a retarded 2-sphere (dark grey) and a rest 2-sphere 
(light grey), defined by an arbitrary timelike worldline in Minkowski spacetime. The 
rest sphere with origin £(to) lies in the orthocomplement of the tangent vector £(to). 
The retarded sphere is the intersection of the rest sphere’s parallel translation along 
£(ro) with the forward lightcone of £(to). The averaging procedure described in the 
text corresponds to sampling the fields of the particle as the radii of these spheres tend 
to zero. Fields on the rest sphere are generated by events on the past history of £(t 2 ). 
By contrast fields on the retarded sphere are all generated by the past history of £(tq). 


zero as well, so this term gives no contribution to the average of the held ()38l) at £(ro). 
Since the second term in (1381) is continuous in a neighborhood of the £(to) its average 
is well defined. Hence the averaged value of (138|) equals the second term on the right 
hand side of (155j) and is also equal to (159j) after an integration by parts. 

An alternative averaging procedure can be defined in terms of rest spheres in the 
orthocomplement of £(to), see Fig. [I] In this case £“(tr(:e)) is not constant on a 
sphere of finite radius r 0 , so the calculation is less convenient; however, in the limit 
of ro tending to zero one finds, again, that the first term on the right-hand side of 
(1381) averages to zero. Thus these two averaging procedures give the same result and 
correspond to surrounding the worldline either by a Bhabha tube or by a Dirac tube. 
See Norton (2009) and Ferris & Gratus (2011) for a similar discussion in the context of 
the standard Maxwell-Lorentz theory with point charges. 
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Motivated by the application of this procedure to inertial motion in example 1 
and the equivalence of the directional-average of (|38j) with (j39j) for general motion, 
we assign to F ab (x) on the worldline the unique value ()39l) for all motions. This is 
intuitively persuasive if one thinks of the point charge as the limiting case of an extended 
charge whose size tends to zero. Directional-averaging is formulated as an axiom in the 
living review on the self-force by Poisson, Pound and Vega, see Section 24.1 in (Poisson 
et al. 2011). 

Given a definition of a finite field F ab (x) at an event x = £(ro) on the worldline, 
the self-force (r 0 ) is defined as the relativistic Lorentz forc^] exerted by F ab (^^(^T o)) on 
the point charge that produces the field: 

/aOo) = qF ab (£(T 0 ))i b (T 0 ). (54) 

The equation of motion for some C° functions £“(r) is assumed to take the form 

m^air) = /f (r) + / a e (r) (55) 


where m denotes the finite inertial mass of the particle and /®(r) is an external 
Minkowski force. If the latter is electromagnetic in origin, /®(r) = q -FJ,(£(t)) £ 6 (t), 
were F^ b solves the Bopp-Podolsky field equations with all sources other than the point 
particle with charge q, which includes the special case of no other sources. If (r) is 
given, (T55li is an integro-differential equation for the worldline £(r). 

In the standard Maxwell-Lorentz theory with point charges, the self-force is infinite; 
therefore, it is necessary to perform a mass renormalization, introducing a “bare mass” 
of the particle that is negative infinite. By contrast, in the Bopp-Podolsky theory there 
is no need or justification for introducing an infinite bare mass. 

From ()39l) the self-force reads 


where 


/aOo) 




d_ 

dC 


w ab (CZ(T 0 )) 


MC/*)dC 

c 


(56) 


VMC,€(to)) = 


(Cb(r 0 ) - ^ b (r))i a (r) - (f a (r 0 ) - £a(r))6,(r) 


Z c (t)(£c( t o) -£ c (t)) 


(57) 


'~=ot(C,£(t 0 )) 

Although (155P contains explicit derivatives of maximal order two, the presence 
of the integral (1561) over the past history of the worldline implies that it cannot be 
solved given only /®(r) and the initial position and velocity of the particle at any 
initial r. Integro-differential equations involving retarded (or memory) effects are not 
uncommon in continuum mechanics and the Maxwell electrodynamics of continuous 
media. Indeed even in vacuo the modifications of the Abraham-Lorentz-Dirac equation 
due to the presence of a background gravitational held yields a similar integral over 
the past history of a charged point particle worldline (DeWitt & Brehme 1960). In 
such problems additional physical criteria motivate analytic or numerical procedures 


* In the appendix, this force is shown to arise from the divergence of the electromagnetic stress-energy- 
momentum tensor associated with the Bopp-Podolsky theory. 
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that can be used to construct solutions. Since in principle the past history of any point 
charge is not empirically accessable, it seems inevitable that similar methodologies will 
be required to determine any future motion uniquely from (1551) given only consistent field 
and particle data on some arbitrary 3-dimensional spacelike hypersurface in Minkowski 
spacetime. These general ideas will be developed elsewhere. 


Example 2: Uniformly accelerated motion 


For a particle in hyperbolic motion with constant acceleration a in its instantaneous 
rest frame, the worldline is given by 


£°(r) = — sinh(ar), ^(t) = - cosh(ar) 
a a 

e(r) = e(r) = 0 . 


In this case, for a point x = £(ro) on the worldline 
2 


reads 


c 2 = 


The self-force ([56]) reduces to 

/aOo) = -q 2 L(r 0 ) 


(cosh(a(r 0 - r)) - 1) . 

50 Mj)d( 


10 


2P,/l + 2y 


(58) 


(59) 


(60) 


which can be expressed in terms of the Bessel functions I\ and K\, 

fa (hr) = - ^2 / i(( a ^) _1 ) | a (r 0 ). 

This result was recently found by Zayats (2014). So the self-force here is manifestly 
finite. 


(61) 


5. Finiteness of the field of a point charge and of the self-force 

In the standard Maxwell-Lorentz theory, the (Lienard-Wiechert) potential of a point 
charge and the corresponding field strength are singular on the worldline of the source. 
By contrast, in the Bopp-Podolsky theory there is a large class of worldlines for which 
the self-force is given by an absolutely converging integral. 

As £(t) is a timelike unit vector, and x — £(r) has Lorentz length (, we may write 

rMC, *)) = cosh x(C, *) <5o + sinhy(C, x) z/ p (C, *) 5 a p (62) 

and 

* a -rMC,*)) = C ^cosh'0(C, x) 6 q + sinh'0(C, x) x) (63) 
where x) and x) are spatial unit vectors, 

Vp{C,x)v p { C,x) = x) /i p (C, *) = 1. 

We introduce the following terminology. 


(64) 
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Figure 5. Worldlinc bounded away from the past light-cone of x (left) and not 
bounded away from the past light-cone of x (right) 


Definition 1. The worldline $, is bounded away from the past light-cone of an event x 
in the future of £ if ?/>(£, x ) stays bounded for ( —> oo. 

Geometrically, £ is not bounded away from the past light-cone of x if and only if 
there is a sequence such that (x — $,(Tk))/(x° — £°(rfc)) approaches a lightlike vector 
for r fc r min , see Fig. [5j 

The notion of being bounded away from the light-cone implicitly refers to a 
particular inertial coordinate system chosen, and it refers to a particular event x. 
However, the notion is actually independent of these choices, as the following proposition 
shows. 

Proposition 1 . The property of the worldline being bounded away from the past light- 
cone of x is preserved if we change the inertial coordinate system by an orthochronous 
Lorentz transformation. If this property is true for one event x in the future of the 
worldline, then the future of the worldline is all of M 4 and the property is true for all 
other events y e M 4 as well. 


Proof. From 


we read that 

x° - £°MC, *)) 


X a - £ a O(C, ®))) [x a - *))) 


= cosh ■?/;(((, a;) 


(65) 


The worldline is bounded away from the light-cone of x if and only if the right-hand 
side is bounded for ( —> oo, i.e., if and only if there exists 5 > 0 such that 

a: 0 — £°(r) 


y- ( xa - £“(r)) (*a - £a(r)) 


< 5 


( 66 ) 
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for r min < r < t 0 with some r 0 . To prove the first part of the proposition, we assume 
that this condition holds in the chosen inertial system. Under a Lorentz transformation, 
x a = A a b x b , the denominator on the left-hand side of ()66l) is unchanged, 

(V - C(r)^j - | a (r)) = (V - f a (T)) (x a - f a (r)) , (67) 

while the numerator changes according to 

x° - £°(r) = A°o(x° - e°(r)) + - ^(r)) . (68) 

As x — £(r) is timelike and future-pointing, 

|x^e(r)|<x 0 -e 0 (r) (69) 


for /i = 1, 2, 3. As a consequence, (1651) implies that 

<k(x° -£°(t) 


(70) 

with some positive constant K. Here we have assumed that the Lorentz transformation 
is orthochronous, A° 0 > 0. From (1661) we find that 

x°-e(r) 


< K 5 


x a — £“(r)) [x a ~ia{r) 


(71) 


for Tm in < t < r 0 which proves that the condition of the worldline being bounded away 
from the light-cone of the chosen event holds in the twiddled coordinate system as well. 
To prove the second part of the proposition, we observe that (163|) implies 


|£-fMC>aO)| 

z° -£°(&(C,x)) 


|tanh^(C,®)| • 


(72) 


Here and in the following, we write |a| = a/ for any a = (a 1 , a 2 , a 3 ). The 
worldline is bounded away from the light-cone of x if and only if the right-hand side of 
(1721) is bounded away from 1 for £ —y oo, i.e., if and only if there is a A with 0 < A < 1 
such that 


|s~ 1*001 

x° — £°(r) 


(73) 


for r min < t < t 0 with some r 0 . Let us assume that this condition holds for some 
particular event x. Let y be any other event and choose a constant y such that 
A < y < 1. Define 

yy° — Xx° — \y — x\ 

t := — - . (74) 

y — A 

Then we have, for all r such that £°(r) < t, 

If-1*001 - y(y° - f°(r)) = \y-x + x- £(r)\ - y(y° - £°(r)) 

< \y-x\ + |x-f(r)| -y(y° -f°(r)) 

< \y-x\ + \(x 0 -e(r))-y(y°-e(r)) 

< \y — x\ + Ax° — yy° + (y — X)t = 0 (75) 
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hence 


-£( r )l 

,0 _ £ 0 ( r ) 


(76) 


for r min < r < fo with some t 0 . This inequality demonstrates that y is in the future 
of the worldline and that the worldline is bounded away from the light-cone of y as 

well. □ 


Because of this result, we may simply say that a worldline is bounded away from 
the past light-cone, without any reference to a specific event x. 

We now show that the field of a point charge is finite if its worldline is bounded 
away from the past light-cone. Using the notation of (1621) and (J63J), the potential (1371) 
reads 

roc (r) a0 + tanhx(C, x) is p ((, x) r) ap ) Ji((/£)d( 

M*) = - 7 / -W- L -t-.(77) 

J 0 cosh-0(C) x) ( 1 - tanh0(C, x) tanhy((, x) yp( C, x) v p { C, x) ) C 


For expressing the field strength tensor at a chosen event x using the notation of (167]) 
and (163]) . we may choose the inertial coordinate system such that £, a (Tn(x)) = Stf. Then 
the electric and magnetic components of the field strength tensor (]38]) read, respectively, 


F 0 a(x) 


2t 2 


+ 


qn a (x) q f°° (tanh^(C, x ) MC, *) ~ tanhx(C, x ) w(C, x )j J 2 (C/Z) d( 


f 2 


tanh0(C, x) tanhx(C, x ) /* p ( C, *) " P ((, x ) C 


,(78) 


Fpa(x) = 


(79) 


Q_ 

e 2 


00 tanh 0 (C,a;)tanhx(C, *) (Ur(C, x)f* P ((, x) - MC, x ) u p((> *)) MC/*■)<% 


1 - tanh'0((', x) tanhy;(C, x ) /^ P (C, x ) u p(Ci x ) ) C 


In a coordinate system with £ a (ro) = the self-force ([56]) is given by 

fa(r 0 ) = (80) 

q 2 r oc ftanh 0 (C,^(r o ))/x ,T (C,^(To))-tanhx(C,^('ro))^ CT (C,C(7'o))) T 2 (C/^)^C 


■ p Vaa 


1 - tanh0(C,^(ro))tanhx(C,^(r o ))/i p (C,^(ro))^(C,l(T-o))) C 

We can now prove the following result. 


Proposition 2. If the worldline £ is bounded away from the past light-cone, the 
integrals on the right-hand sides of (1771) . (1781) and (|79l) are absolutely convergent for 
all x E M 4 ; the integral on the right-hand side of (180]) is absolutely convergent for all 

Tmin G M) <: t 7" max . 


Proof. As the integrands in (177]) . (178]) . (179]) and (180]) stay finite for ( —> 0, we only have 
to verify that they fall off sufficiently quickly for ( —y 00 . We first observe that tanh a 
increases from —1 to 1 if a varies from — 00 to 00 , and that |/TV p | < 1. Therefore, 
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the condition of ip being bounded implies that tanh ip tanh y /i p // p is bounded away from 
1. Thus, in each of the four equations (177]) . (178]) . (179]) and (180]) the modulus of the 
integrand is bounded by a term of the form K\Jk((/£){/( where K is independent of 
C and k is either 1 or 2. As |Jfc(C/^)| falls off like C -1 ^ 2 for C ~t °°> this guarantees 
absolute convergence of the integral. □ 

This proposition demonstrates that, in particular, the self-force is finite for a large 
class of worldlines. Actually, the requirement of the worldline being bounded away 
from the past light-cone is sufficient but not necessary for finiteness of the self-force. 
The following proposition shows that there is another class of worldlines, including ones 
which are not bounded away from the past light-cone, for which the self-force is finite. 

Proposition 3. Assume that the worldline £ is confined to a two-dimensional timelike 
plane P in Minkowski spacetime. Then the integral on the right-hand side of (]80l) is 
absolutely convergent for all r E ain < To < T max . 

Proof. As the worldline is in a timelike plane, the unit vectors £( r o)) an d 

/hr(C,£(7b)) must be constant and equal or opposite, u a = ±/i CT . Then (180]) simplifies to 

/a(ffi) = tanh^(C, g( 7 o)) T x(C, €(t 0 ))) ^ . (81) 

As |tanh(a:)| < 1, the modulus of the integrand in (]8H) is bounded by |^(C/^)I/C- As 
| J 2 (C/^)| falls off like C _1/2 for C t oo, this guarantees absolute convergence of the 
integral. □ 


Example 3: A worldline with diverging self-force integral 


From Propositions [2] and [3] we know that the self-force is finite if the particle’s worldline 
is bounded away from the light-cone or if it is contained in a timelike plane. Actually, 
the proof of Proposition [3] can be generalized to the case that the worldline, rather than 
being confined to P, approaches P sufficiently quickly for r —$■ r min . This leaves only 
a class of rather contrived motions for which the self-force integral could diverge: The 
worldline must approach the light-cone for r —> r m ; n with a sufficiently large tangential 
velocity component. In this section we present such an example for which the self-force 
integral, indeed, diverges at one point. 

We find it convenient to give the worldline in terms of a past-oriented curve 
parameter 7 which is not proper time, 

£ a (l) = b{ - V /l + s (7)<5o + \/^( 7 ) ( cos (/^( 7 ))5“ + sin(/ 3 ( 7 ))h“)} (82) 

where 


and 

£(7) 


s(7) 


7 

i tanh ( 7 /€) 


vTrC' J F(v / 2s(7)/7r) - y/7rS F (yj2s('y)/n) - 


sin(s( 7 )) + cos(s( 7 )) 

V 2s (V 


(83) 


(84) 











Figure 6. The worldline (15211 approaches the past light-cone (dotted circle) with an 
oscillatory tangential velocity component 


(85) 


Here Cf and Sf are the Fresnel-C and Fresnel-S functions. Note that 

s'( 7 )sin(s( 7 ) + tt/4) 

H,) = —17733—• 

£ is an analytic timelike curve that approaches the past light-cone with an oscillatory 
tangential velocity component, see Fig. El The curve is, indeed, everywhere timelike as 
can be seen from 

^“(7)^ 6 (7) , , 7 2 s'(7) 2 


Vaf 


dy dy 


= - 1 + 


< -i + i + 7 = 

4 16 


4s( 7 )(l + s( 7 )) 16 s( 7 ) 2 

11 


+ sin 2 (s( 7 ) + 7 t/4 ) 


16 


( 86 ) 


It can be shown that the 4-acceleration of $, is bounded and that its future is all of 
Minkowski spacetime. 

We show that for this worldline $, the electromagnetic field (1381) diverges on a 
hyperplane implying that the self-force becomes infinite at one instant. To that end we 
have to rewrite the integral in (l38|) as an integral over y and to investigate the behavior 
of the integrand for 7 —> 00 . We first observe that, if 7 tends to 00 , 


5(7) =j + 0(y~ n ), s'(y) = j + 0(y ~ n ), 


/%) = 


VI 


Sill 


+ 


7 


3/2 


+ 0(y~ n ) 


(87) 


for all n6N. Moreover, from the standard asymptotic formulas for the Fresnel functions 
we find 

cos ( j 3 ( y )) = 1 + 0 ( 7 -3 ), 8 m ( j 3 ( y )) = 0 ( y ~ 3/2 ). (88) 

With the help of these formulas we find that the parameter ( which is used as the 
integration variable in (l38]i is related to our curve parameter 7 by 

3/2 

< 2 = -(/- f (7 ))d+ - u-r)) = 7 2 + 2-b=-(x° + x 1 ) + 0 ( 7 "), (89) 
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hence 


c = T+ y|( I o +lI )_(^±f2! + 0(T - I / 2 ) 


(90) 


Inserting this expression into the well-known asymptotic formula for the Bessel function 
h yields 


J?(-) = 


IHK)+°(c- 3/2 ) 


(91) 


2<. (1 , v ^j{x a + x y ) (xo + njy 

- Sm ( ~i) + — H-7T/5-75- ) + 0(7 J 

7T7 V £ 4 £ 3 / 2 £2 / 


After these preparations, we are ready to evaluate the integral in (1381) . If we use 7 as the 
integration variable, writing 7 r(x) for the parameter value that corresponds to tr(x) 
and thus to C = 0 , this integral reads 

- (to - 6(7))“ - (*. - 6,(7))“) j 2 ( C/ q dc 


hb 


'7r(*) 


- « 7 » C <*7 

- (fo - &(7))*“ - (*. - &(7))“ i ) j 2 ( C/ q ,i C 


'ir( x ) 


1r(, x ) 


^ c?7 


C ^7 


d7 (92) 

d'y 


(x b - 6 ( 7 )) 


^a(7) 

d'y 


[x, 


a sa 


,, 22 rfa(7)A ^(c/«), 

iM) —) 


We evaluate this equation for a = 2 and 6 = 0. As 


Oo - £ 0 ( 7 )) 


2(7) / * / 77^0(7) -^ 3/2 


(®2 -6(7))- 


-7 


c?7 v ~^ ^ w " d7 4v^ 

we find with our asymptotic formula for ^(C/^) from above 


.sm(l + 0 + O( 7 '7. 


(93) 


ho — 


:sm( ] + f 


h r( x ) 


2a/27t 


Sill 


7 


G + j + ^)-^) + ot^)U 


4^2^ 


'7h(*) 


cos 


^(^T^ 1 ) (x° + x 1 ) 2 


-cos 


£ 3/2 £2 

27 7r 77^° A^ 1 ) (x 0 -^ 1 ) 2 


(94) 


d 7 

— + ... 

7 


£ 2 £ 3 / 2 i 2 

Here and in the following, the ellipses indicate a term that is finite for all x. The integral 
over the second term is finite for all x. If we decompose the remaining integral into an 
integration from jn(x) to l and an integration from t to infinity we find 
'h~f(x° + x 1 ) (x° + x 1 ) 2 \d'y 


ho — 


4\/27r 


cos 


A\/27t { 


I- 


2 cos 


£ 3/2 

(x° + X 1 ) 2 


£ 2 


7 


£2 


Ci 


x 


• Sill 


(x° + a; 1 ) 2 
£2 


+ x 1 ! 3 / 2 

~W 2 

7T — 2 Si 


(95) 


\x° + X 1 ! 3 / 2 

W 2 
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Here Ci and Si denote the cosine integral and the sine integral, respectively. As the 
cosine integral diverges logarithmically if its argument approaches zero, we have found 
that the (20)-component of the electromagnetic held according to ()38l) is given by 



(96) 


i.e., that this held component diverges on the lightlike hyperplane x 0 + ad = 0. This may 
be interpreted as a shock front propagating at the speed of light. The same divergence 
is found, by a completely analogous calculation, for the component F‘2l — —Fi2, 



(97) 


while all other components are finite everywhere. As a consequence, the self-force 


becomes infinite at the instant when the charged particle crosses the hypersurface 
x° + x l = 0 which happens at the origin of the coordinate system. This means that at 


this instant an infinite external Minkowski force is necessary to keep the particle on its 
prescribed worldline. 

Note that the divergence is logarithmic in the neighborhood of a lightlike hyperplane 
and thus rather mild, since any timelike worldline crosses it at most once. It is true 
that a charged particle would experience an infinite relativistic Lorentz force at the 
instant when it crosses the hypersurface x° + x 1 =0. However, as f y log|x| dx is finite¬ 
valued and continuous at y — 0, the particle’s velocity would still be hnite-valued and 
continuous, i.e., the particle’s worldline would still be a C l curve. 

6. The Abraham-Lorentz-Dirac limit 

In the standard Maxwell-Lorentz theory with point charges, i.e., for i —)■ 0, the self-force 
becomes infinite in ()54lh Dirac’s solution to give a meaning to the equation of motion 
(T55|) in this case was to assume that the inertial mass became negative infinite in order to 
cancel the infinite contribution from the self-force. After this cancellation, one ends up 
with the Abraham-Lorentz-Dirac equation which involves a renormalized (or dressed) 
mass which is positive and hnite. It is interesting to see how the Abraham-Lorentz-Dirac 
equation is reproduced from the Bopp-Podolsky theory in the limit 7 —>■ 0. To that end 
we substitute in (|56|) the integration variable ( = £a, 



(98) 


where 



(99) 


With 


x(0) = i, :v'(0j - -1 


( 100 ) 
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two times integrating by parts yields 

fa ( A)) = - |^la(To) + ( L(a) + L(t 0 ) ^(tq) £ b (t 0 )) (101) 

/ °° ffiW 

-^(£a^(r 0 )) X (a)da. 

The first term diverges for i —» 0. Following Dirac’s idea of mass renormalization, the 
parameter m must depend on i and become negative infinite such that the “dressed 
mass” 

2 

™ = I™ (m(£) + (102) 

remains hnite and positive. In this limit, the equation of motion reads 

fnia(r) = la (t) + L(t)£\t) l b (t)) + /®(t) (103) 

( roc \ 

V£ 6 (a) —^y(l<J,£(TQ))x(v)daj . 

If the integral is bounded, the last term vanishes and we get the Abraham-Lorentz- 
Dirac equation. From (jU9|) we find, with the help of the well-known asymptotic formula 
for the Bessel function Ji, that y(cr) = 0{a~ 3 ^ 2 ) for a —> oo. So the integral in (I103jl is 
certainly bounded if d 3 W a b/d( 3 is bounded for ( —> oo. A sufficient (but not necessary) 
condition is that £ is bounded away from the light-cone and that all components £ a , £“, 
£“, and £“ are bounded. 

7. Conclusions 

In this paper we have demonstrated that in the Bopp-Podolsky theory a self-force of a 
charged point particle can be defined by an integral that is absolutely convergent for a 
large class of worldlines on Minkowski spacetime. We have also provided a (contrived) 
example where its electromagnetic field diverges on a lightlike plane, so the self-force 
diverges at one point of the worldline. However, even in this case the electromagnetic 
field is a locally integrable function (i.e., a regular distribution), yielding a C 1 solution 
to the equation of motion (J55]). This is to be contrasted with the standard Maxwell- 
Lorentz theory in vacuo where the self-held of a charged point particle is infinite at 
every point of the particle’s worldline, and the singularity is so bad that the held energy 
in an arbitrarily small ball around the charge is infinite necessitating classical mass 
renormalization. In the Bopp-Podolsky theory there is no need for mass renormalization; 
the equation of motion (1551) is an integro-differential equation making sense with a hnite 
inertial mass m. 

It should also be emphasized that in the Born-Infeld theory, which is to be viewed as 
a natural rival to the Bopp-Podolsky theory, virtually nothing is known about hniteness 
of the self-force of an accelerated particle whose worldline may approach the light-cone. 
So it seems fair to say that, at least in view of the motion of charged point particles, 
the Bopp-Podolsky theory is in a more promising state. 
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Some questions remain open. It is clearly important to establish a strategy for 
solving (1551) . This may include formulating a consistent Cauchy problem that may 
include either background or dynamic fields satisfying the Bopp-Podolsky held equations 
coupled to the particle. Furthermore, it would be desirable to have a proof that all 
worldlines solving ()55l) are at least C 1 . It is also important to demonstrate that the 
equation of motion (T55l) with vanishing external Minkowski force is free of run-away 
solutions. Partial results in this direction have been found by Frenkel & Santos (1999), 
but a general proof is still missing. If these questions can be satisfactorily resolved, the 
Bopp-Podolsky theory offers a physically acceptable framework in which to explore the 
classical electromagnetic back-reaction on charged point particles. 
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Appendix 

The formulation of a stress-energy-momentum tensor for the theory discussed in this 
paper appears to have had a chequered history. Bopp (1940) writes down a stress-energy- 
momentum tensor that is obviously based on the decomposition (ITU of the potential, 
but no derivation is given. Podolsky (1942) also writes down a stress-energy-momentum 
tensor with a promise to derive it in (Podolsky & Kikuchi 1944) from arguments based 
on a canonical approach. In our view this did not succeed and the further derivation 
in (Podolsky & Schwed 1948) lacked transparency. To our knowledge there has been 
no subsequent attempt to derive any stress-energy-momentum tensor appropriate to 
the theory under discussion. In view of these comments it may be of value to put the 
matter into a modern perspective by offering a derivation based on metric variations of 
the Bopp-Podolsky action. This requires formulating the theory on a curved spacetime 
manifold. 

The natural tools for this purpose exploit the exterior calculus of differential forms 
using properties of the Hodge map * associated with the spacetime metric g and the 
nilpotency of the exterior derivative d. For background material on exterior calculus 
we refer to Straumann (1984) whose sign and factor conventions we adopt. The Bopp- 


(Aem — A A J) 
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Podolsky action reads 



(104) 


with 



(105) 


where F = dA and G = d * F are smooth. Here we assume that J is a prescribed 
(non-dynamical) smooth current 3-form that satisfies the conservation law dJ = 0 on 
A4. In the body of the paper we have restricted ourselves to the case of a flat metric 


and used inertial coordinates. Then the action (1104)1 reduces to (J3J) where the current 


4-vector j = j a d a is related to the current 3-form by g(j, ■) = *J. 

A direct route to the symmetric dynamical (Hilbert) stress-energy-momentum 
tensor, associated with Aem, is obtained by making compact variations of the metric 
tensor in Aem- Such variations can be induced by making independent variations e° in 
a local g— orthonormal coframe {e a }, a = 0,1, 2, 3 since in such a basis g = r] a be a <£> e b . 
Such variations give rise to a set of 3-forms r a defined by 



(106) 


and a stress-energy-momentum tensor T = T a i,e a <£> e b with components 
Tab = Vbc * (r a A e c ). 

The covariant divergence of T then follows as 
V • T = E'DrJe- 


(107) 


(108) 


where D denotes the covariant exterior derivative, see e.g. Benn & Tucker (1988). Since 
T is symmetric 


Dr a = dr a - i a de b A n . 


(109) 


If we make compact variations of the potential A in S, rather than of the metric, 


we obtain the field equation of the Bopp-Podolsky theory. We can derive the r a and 


the field equation in one go if we allow for partial variations of the potential and of the 
metric simultaneously. (Note that the current J is assumed to be given and is kept fixed 


during the variation.) Then the total variation of the Lagrangian 4-form is written 



( 110 ) 


87r 87T 


For calculating the extremum of A we use the formula (Dercli et al. 2007) 


(*T)'= e' c A i c (*T) - *(e’ c A i c T) + 


( 111 ) 


that holds for any p —form T, where i c denotes the contraction operator (or interior 
derivative) with respect to the vector field X c defined by e a (X c ) = <5“. Moreover, we use 
standard rules of exterior calculus, such as T A = $ A for any p— forms T, $, the 
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graded derivative property and nilpotency of d, and the commutativity of d with the 
variations. Thus 


(F A -kF)’— 2d(A’/\*F) + 2A’Ad*F + e’ c /\ (F A i c * F — i c F A F)(112) 


and 


{G A *G)' = 2 d ( (e' c A i c -k F — *(e c A i c F) + A ★<?) 

— 2 A A d -k d -k G — e r A (2i c -kFAd~kG — 2 i c F A -kd -k G 


+ G Ai c *G + i c G A *G 
Inserting (11121) and (11131) into (11101) yields 


where 


A' = + e' c A r c + A' A —dH — A'AJ 

An 

An<& = A’ A H — t 2 i^e' c A i c * F — -k(e’ c A i c F ) + A -kG , 

H = -kF -\~^~kd-kd-kF 


(113) 

(114) 

(115) 

(116) 


and 

87 tt c = F Ai c ~kF — i c F A-kF(^G Ai c ~kG-\~icG A-kG — 2i c F A-kd-kG-\~2i c -kF Ad-kG^j .(117) 

The field equations are determined by requiring that the action is stationary for 
partial variations of the potential only (i.e., e c = 0 ) that are compactly supported (i.e 

dH = 4wJ. (118) 

Similarly for partial variations with A = 0, the r c from (11171) give, via (1107(1 . 
the dynamical stress-energy-momentum tensor of the Bopp-Podolsky theory which is 
automatically symmetric. A lengthy but routine calculation of Dr a then shows, with 
the aid of the field equations (II 18(1 . that the divergence (11081) yields the relativistic 
Lorentz force, 

V-T = F{-,j). (119) 

A less lengthy approach to derive (11191) can be based on the use of a one-parameter 
family of diffeomorphisms on the spacetime domain M. generated by any compactly 
supported vector field X. Then in terms of the Lie derivative Lx 

Lx Aem = dA>x + Lx& c A r c + -—Lx A A dH (120) 

47T 

with the same r c and the same H as in (j 11 4(1 and 

An ( L y = L x A AH -f(^L x e c Ai c *F -*(L x e c Ai c F)+*L x F^ A*G.(121) 
For the derivation of (11201) we have used the identity 


L x {kAi) — L x e c A i c (*T) — -k{L x e c A i c T) + kL x A> 


( 122 ) 






















which holds for any p-forrn T on M. With the relations Lx = ixd + dix and <7 Aem = 0 
(since A EM is a 4-form on spacetime), (11201) yields the identity 

doL X = fix (123) 


where 


and 


olx = ^xAem — ixe c A r c — ixA A —dH — <3>x 

47T 


fix — i\de c A t c — ixe c A dr c + i\F A —dH . 

47r 


(124) 


(125) 


Integrating (11 23 ft over an open domain U containing the support of X and using Stokes’s 
theorem yields 


fix — da x — ol x — 0 


<u 


<u 


>dU 


since ax\au — 0- Furthermore since fix has the linearity property fifx = f fix for 
any smooth function /, it follow^ that fix = 0 (Gratus et al. 2012). Finally choosing 
X = X a and using the field equation (11181) together with (11091) gives 


Dr a = F(X a , •) A J (126) 

which is equivalent to (11191) . using (I108jl . From this derivation one concludes that, for 
any diffeomorphism and gauge invariant action 

S[A,J,g\= I A (A,J,g) 

JM 

constructed from 

A(A, J, g) = A f (F, g) — AA J (127) 

r 

and F = dA equation (11261) is satisfied when J and g are background fields with A I (F, g ) 
arbitrary. However for consistency this requires that J be a prescribed exact 3-form. 
On a topologically trivial spacetime domain A4 this is implied by dJ = 0. 

If the background metric is flat and inertial coordinates are used, (11181) reduces to 
(151); in this case the components of the stress-energy-momentum tensor associated with 
(11171) are 

4rrT cd = -J ab F ah r }cd - F db F c b + f (d b F bc d a F ad - l -d a F ab d,F%J) (128) 
- e ( F cb d a d a F b d + F db d b d a F ac + F cb d b d a F ad + F ab d e d e F abVcd ^j 

where we have used dF = 0. In this flat metric with the field equations (J5J) , the stress- 
energy-momentum tensor (11281) satisfies d c T cd = F db j b which is the flat-space coordinate 
version of (11191) . The stress-energy-momentum tensor (I128jl coincides with that written 
down by Podolsky (1942) and subsequently quoted by Zayats (2014). The derivation 
above shows that, in general, it coincides with the Hilbert stress-energy-momentum 
tensor derived from metric variations of the action (11041) . 


jt Note however since ajx f olx for arbitrary /, one cannot similarly conclude that ax = 0. 
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Finally we point out how different definitions of the stress-energy-momentum tensor 
for the Bopp-Podolsky theory are responsible for the historic problems outlined in the 
beginning of this appendix. In general (J9]) becomes 

A = -f*G, A = A- f -kG. (129) 

and the Lagrangian (1104j) can be rewritten as 

G 

A = A 1 + — d(k d-k F A kF ) (130) 

47r 

where 

A 1 = -5 -gL 4 A kdA - (dA A kdA + \-A A kA J - (A - A) A J . (131) 

87r 87T y G J 

As A and A 1 differ only by an exact form, they lead to the same field equations and to 
the same Hilbert stress-energy-momentum tensor. Indeed, varying the action 


S'[A As] 



j ( — dA A 'kdA — — ( dA A ~kdA H—— A A krA 

Jm V8tt 8tt V P 


(132) 

(A-A) A J ) 


with respect to A and A respectively yields the field equations 

dkdA = 4:7rJ 1 d k dA — — k A = AnJ. (133) 

These equations imply that A necessarily satisfies the Lorenz gauge condition, dkA = 0. 

The field equations (1 133p reduce to (J7J) and (jHJ) in flat spacetime using inertial 
coordinates and the Lorenz gauge condition on A. 

The g —orthonormal co-frame variation of A 1 gives rise to the Hilbert stress forms 

8ttt c = dAAi c kdA—i c dAAkdA—dAAi c kdA+i c dA/\kdA+ — (^AAi c kA+i c AAkA' S j .(134) 

Substituting from (1129ft . we see that these stress forms indeed coincide with (11171) . On 
flat spacetime in inertial coordinates, (1134ft yields the stress-energy-momentum tensor 

47rT cd = A ab F ah ri cd - F ac F a d - A ab F ab r, cd + F ac F a d + 1 (A c A d - A a A a rj cd ) (135) 

where F a b = d a A b — d b A a and F ab = d a A b — d b A a . This is the stress-energy-momentum 
tensor given by Bopp (1940), cf. again Zayats (2014). 

Furthermore the same stress forms (11341) arise using the canonical Belinfante- 
Rosenfeld symmetrization of the Noether current associated with the Lagrangian (1131ft . 
By contrast the canonical Belinfante-Rosenfeld procedure applied to the action ([2]) does 
not yield the stress-energy-momentum forms (11341) . This illustrates the fact that two 
Lagrangians that differ by an exact form, while yielding the same Hilbert stress-energy- 
momentum tensors, in general, yield different stress-energy-momentum tensors following 
the canonical Belinfante-Rosenfeld procedure. 
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